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ABSTRACT 

An  investigation  is  made  of  the  wave  forecasting  method  of  Pierson, 
Neumann  and  James  of  New  York  University,  published  as  II,  0.  Pub.  No,  603. 
The  oooan  "swell"  speotral  energy,  as  predicted  by, the  application  of  this 
forecast  method  is  compared  to  the  observed  spectral  enorgy  present  as 
estimated  from  ooean  wave  recording  devices  installed  at  the  surface  and 
at  a  depth  of  80  feet,  in  the  coastal  waters  off  Davenport,  California. 

The  comparison  of  the  observed  and  forecast  spectral  energy  for  nine 
verification  times  taken  at  twelve  hour  intervals  from  28  January  1954  to 
1  February  1954  suggested  the  following  conclusions: 

1.  For  the  conditions  of  the  invest! gat ion#  it  was  the  opinion  of  the 
authors  that  the  "swell"  forecast  was  operationally  useful.         ** 

2.  The  data  of  this  investigation  suggested  the  possibility  that  the 
foreoast  energy  spectrum  is  too  high  in  the  low  and  middle  frequencies  and 
too  low  in  the  high  frequency  bands  for  high  wind  speeds. 

3«  The  data  of  this  investigation  indicated  that  at  certain  verifica- 
tion times,  energy  reductions  due  to  cross-sea  interference  were  important. 

The  authors  wish  to  express  their  appreciation  for  the  advice  and 
enc ouragement  of  Professor  J.  B.  Wickham  under. v/hose  direction  this  study 
was  carried  out. 

Appreciation 'is  also  expressed  to  Professor  R.  L.  Wiegel,  of  the 
Institute  of  Engineering  Research,  Berkeley,  California,  who  made 
available  the  wave  records  upon  whioh  this  study  is  based. 
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CHAPTER  1 
INTRODUCTION 

The  value  of  reliable  techniques  for  forecasting  ooean  wave  proper- 
ties for  seagoing  and  coastal  operations,  particularly  military  operations, 
is  obvious* 

This  problem  has  engaged  the  interest  and  attention  of  some  of  the 
foremost  scientists  of  various  fields  over  the  years.  Some  of  the  contri- 
butions have  been  notable. 

The  theories  of  classical  hydrodynamics  were  applied  by  Lamb  QlJ 
in  1879,  which  yielded  important  properties  of  oertain  relatively  simple 
sine-wave  forms. 

In  1933,  CornishQfJ,  published  the  results  of  a  large  number  of 
carefully  made  visual  observations  of  wind-waves,  both  in  deep  water  and 
near  the  shore,  which  were  the  first  of  that  kind  and  which  undoubtedly 
stimulated  further  work  of  this  nature. 

The  development  of  continuous  wave  recording  devices  [_Z J  expanded  the 
wave  data  available  and  made  possible  the  many  statistical  studies  which 
have  been  carried  out. 

Sverdrup  and  I.Iunk  QQ   published  a  wave  forecasting  method  in  1947, 
for  predicting  the  generation,  decay,  and  transformation  of  the 
"significant"  wind-waves. 

Bret Schneider,  in  1951,  published  a  modification  forecast  [V]  of 
the  Sverdrup  and  T.lunk  technique  based  on  more  observational  data. 

In  1953,  Pierson,  Neumann  and  James  of  New  York  University  pub- 
lished a  new  theory  and  corresponding  foreoast  technique  (1Q  •  This 
theory  and  method  are  in  many  respects  unique  and  depend  largely  on  some 
assumptions  which  have  had  limited  testing. 

1 


Sinoe  this  last  theory  and  technique  has  been  adopted  for  use  by  the 
United  States  Naval  Service,  a  study  and  evaluation  of  the  method  is  the 
purpose  of  this  paper,  A  study  of  this  nature  presents  many  difficulties 
and  requires  data  and  facilities  available  at  very  few  inst itut ions •  These 
facts  explain  the  paucity  of  published  results  of  such  studies*  To  the 
knowledge  of  the  authors,  no  such  studies  and  evaluations  havo  been  made 
previously  of  the  Pierson,  Neumann  and  James  wave  forecast  method,  hence- 
forth rof erred  to  as  the  N.  Y.  U.  method,  for  Pacific  Ocean  swell,  although 
the  waves  forecast  in  a  generating  area  have  been  studied  [l6j. 

The  study  divides  itself  naturally  into  four  distinct  phases: 
1*  A  wave  forecast  for  a  given  period  of  time  at  a  specified  location, 
using  the  11.  Y.  U.  method, 

2.  The  analysis  of  wave  records  from  a  wave  recording  device  for  the 
same  location  and  same  time  to  obtain  the  actual  wave  conditions. 

3.  The  comparison  of  the  forecast  and  observed  wave  properties. 

4.  The  results  or  conclusions  to  be  drawn  from  the  comparison. 

The  data  and  facilities  necessary  to  carry  out  this  study  are  listed 
below  with  their  sources: 

1.  A  forecast  looation  at  which  suitable  wave  records  have  been 
kept  for  a  convenient  period  of  time.  The  site  chosen  was  Davenport, 
California,  where  the  University  of  California,  Berkeley,  had  a  series  of 
recording  devices  [7 3  installed  from  the  period  of  Nov.  1952  to  Dec.  1954. 
Wave  records  from  Jan.  1954  to  Feb.  1954  were  made  available  to  the 
authors  for  use  during  this  study. 

2.  A  rapid  computing  device  for  the  analysis  of  the  wave  records  [V] 
to  obtain  the  distribution  of  wave  energy  as.  a  function  of  frequency  Q?j. 
The  U.  S.  Naval  Postgraduate  School  is  equipped  with  a  CRC  102-A  medium- 
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speed  general-purpose  electronic  computer  which  was  available  for  this  purpose, 

3.  Synoptic  surface  weather  maps  for  the  period  of  and  proceding 
the  forecast  time  for  the  purpose  of  preparing  the  detailed  wave  foreoast* 
Such  maps  were  made  available  by  the  Department  of  Mechanical  Engineering 
Yfave  Research  Laboratory,  University  of  California,  Berkeley. 

4,  A  submarine  contour  ohart  for  the  area  adjacent  to  the  forecast 
location.  U.  S.  Coast  and  Geodetic-  Original  Survey  Sheet  Ho»  5266  of 

25  April  1935,  was  made  available  by  the  Hopkins  Uarine  Station,  Pacific 
Grove,  California. 


CHAPTER  2 
TIIE  THEORY  AND  THE  METHOD 

The  objective  of  this  chapter  is  to  discuss  the  portions  of  the 
N.  Y.  U.  Method  £6j  and  theory  -which  relate  diroctly  to  this  investiga- 
tion. 
2A       The  Spectrum  of  Wind-Generated  Waves 

The  primary  basis  for  this  forecast  method  is  the  assumed  distribu- 
tion of  ocean  wave  energy  for  a  given  wind  speed  as  a  function  of  wave 
frequency  in  a  wind-wave  generation  area,   or  the  energy  spectrum.     This 
distribution  is  represented  mathematically  by: 

g2 
CA(*T  «    ,     °       a  e  "  27r2f  2V2  U) 


T^F?&       27T2f2v2 


Where  o  is  a  constant  of  proportionality  and  the  other  terms  are  defined 
on  page  VI.  Figure  1  shows  a  plot  of  this  function  for  wind  speeds  of  20, 
30,  and  40  knots,  and  tabulated  values  of  [A(f  )J   as  a  function  of  wind 
speed  and  frequency  are  given  in  Appendix  I. 

The  form  of  this  function  was  assumed  by  Ileumann  (jLOj  based  on 
theoretical  considerations  and  empirical  evidence.  To  a  large  measure, 
the  accuracy  of  this  method  will  be  determined  by  how  closely  the  assumed 
spectrum  approximates  the  true  spectrum.  Such  determination  will  be 
possible  only  as  increased  test  results  become  available. 
2A1  Energy  in  the  Spectrum 

The  integral  of  Equation  (l)  for  a  given  wind  velocity  with  respect 
to  frequency  is  proportional  to  the  energy  of  the  included  spectral  com- 
ponents. This  measure  of  spectral  energy,  between  the  frequencies  f 
and  fa  $   *s  defined  as  E  and  is  given  by 


■  ft  «>] 


df  (2) 


fl 


and  has  the  dimensions,  length  squared. 

Well-tested  statistical  relationships  [lljl  between  E  and  wave  heights 
have  been  demonstrated.  These  relationships  are  given  in  Table  1  and  are 
used  as  a  basis  for  all  wave  height  determinations,  E  being  the  forecasted 
parameter. 
2A2  Properties  of  the  Spectrum 

The  wind-generated  sea  for  a  given  wind  speed  is  said  to  be  fully 
developed  when  all  possible  wave  components  in  the  spectrum  between  f  a  0 
and  f  =  00  are  present  with  their  maximum  amounts  of  spectral  energy. 
This  state  is  reached  only  after  the  wind  has  blown  for  a  sufficient 
length  of  time  ("duration")  and  over  a  sufficiently  long  distance  ("fetch"). 

If  the  wind  duration  or  fetch  is  insufficient  for  the  fully  developed 
state,  a  lower  energy  level  will  result.  The  speotral  components  form  in 
the  high  frequency  end  of  the  spectrum  initially  and  the  lower  frequency 
components  are  added  progressively  and  grow  as  the  wind  duration  increases 
until  full  development  occurs,  providing  the  fetch  is  not  limiting  and  the 
wind  continues.  One  important  feature  of  this  theory  shows  the  partially 
developed  spectra  for  limiting  conditions  of  fetch  and  duration  for 
various  wind  speeds. 
2A3  The  Co-oumulative  Spectrum 

The  energy  spectrum  is  oonvorted  to  a  more  convenient  form  by 
integrating  it  over  the  entire  frequency  range.  This  procedure  is  illus- 
trated in  Figure  2.   It  will  be  noted  that  the  ordinate  of  the  cumulative 


integral  can  be  read  diroctly  in  values  of  E,  which  is  proportional  to  all 
of  the  wave  energy  in  spectral  components  equal  to  or  groater  than  the 
abscissa,  f.  The  functional  relationships  between  energy,  wind,  duration 
and  fetch  may  also  be  included  on  the  cumulative  spectrum* 

Figure  3  and  Figure  4  illustrate  the  IT,  Y,  U,  co-cumulative  spectra 
given  as  functions  of  duration  and  fetch  respectively  for  various  wind 
speeds.  These  are  the  curves  normally  used  in  the  N,  Y,  U,  method, 
2B   The  Wave  Forecast  in  the  Generating  Aroa 

The  foregoing  are  required  for  forecasting  the  wavo  conditions  in  a 
wind-wave  generating  area.  Knowing  the  wind  speed,  the  wind  duration,  and 
the  fetch  length,  the  E-value  and  the  range  of  wave  periods  with  import- 
ant energy  may  be  obtained  diroctly  from  the  co-cumulative  spectrum.  The 
sea  in  the  generating  aroa  has  an  irregular,  short-crested,  appearance  as 
the  spectral  components  randomly  amplify  and  damp  each  other.  The 
spectral  components  also  move  in  a  variety  of  directions,  as  the  wind 
action  is  not  uniform,  but  of  varying  strength  and  direction, 
2C   The  "Swell"  Forecast 

When  wind- gene rated  waves  leave  the  generating  area,  a  transforma- 
tion takes  place  in  the  appearance  of  the  sea  surface  as  the  wind  no 
longer  acts  on  the  wave  form.  As  the  various  spectral  componont  wave 
trains  leave  the  wind  area,  the  apparent  waves  becomo  regular  and  long 
crested,  and  are  termed  "swell."  When  "swell"  arrives  at  points  distant 
from  the  generating  area,  the  waves  have  only  a  fraction  of  the  energy 
originally  present  in  the  fetch, 
2C1  Assumed  Processes  Involved  in  Energy  Reductions 

Pier son,  Neumann  and  James  account  for  energy  reductions  at  points 
distant  from  the  fetch  by  two  processes,  namely 
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1.  Dispersion 

2.  Angular  Sproading 
202     Dispersion 

When  the   speotral  components  leave  a  generating  area,  the  associated 
energy  front  of  each  travels  at  different  group  velocity*     The  group 
velocity  for  any  spectral  frequency  is  given  by 

G.    V.   a   1.515  T  (3) 

Dispersion,  then,   is  defined  to  be  "the  spreading-out  effect  caused  by  the 
different  group  velocities  of  the  spectral  frequencies   in  the  original 
disturbance  at  the   source."   [jfj    Therefore,  to  determine  for  a  given  time 
the  total  spectral  energy  at  a  specified  distance  from  the  fetch,  based  on 
the  effects  of  dispersion,    it   is  necessary  to  compute  the  range   of  fre- 
quencies which  could  be  present  at  that  time.     This  may  be  accomplished  by 
considering  the  distance  the  energy  must  travel,  the  group  velocity,   and 
the  time  the  particular  spectral  component  was  formed  as  a  function  of 
wind  speed,   duration,   and  fetoh. 

The  N.  Y.  U.  method  includes   some  basic   "filters"  which  are  essen- 
tially pre-arranged  formulations  of  these  variables,   for  certain  recurring 
situations.     The  use  of  these  "filters"  permits  the  calculation  of  a  Ty, 
upper  period,  and  Tt,   lower  period.     Ty  and  T^  define  the  upper  and  lower 
limits  of  the  spectral  energy  present  at  the  forecast  time  and  location, 
for  a  given  wind  speed. 
2C3     Angular  Spreading 

The  considerations  of  paragraph  2C2  do  not  take  into  account  the 
effect  of  variability  in  the  direction  of  wave  motion  and  the  fact  that 
the  waves  are  short-crested.     The  component  wave  trains  vary  not  only  in 
frequency  but  also  in  direction  within  the  area  of  generation.     Thus,   as 
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the  spectral  components  leave  the  generating  area,  they  spread  out 
angularly  as  -well  as  "disperse." 

The  assumed  distribution  of  relative  wave  energy  as  a  function  of 
diroction  is  illustrated  in  the  lower  curve  of  Figure  5,  Since  the  area 
under  this  curve  is  unity,  when  integrated  cumulatively  over  direction  it 
gives  that  fraction  of  the  energy  to  be  found  between  any  two  directions 
(measured  from  that  of  the  mean  wind).  The  ordinate  of  the  upper  curve  in 
Figure  5  is  then  expressable  in  percent.  Dispersion  and  angular  spreading 
computations  are  illustrated  in  chapter  3,"~ 

Figure  6  illustrates  the  measurement  of  the  angles  ©^  and  Q§,   which 
are  the  parameters  used  in  conjunction  with  the  upper  curve  of  Figure  5 
for  evaluating  the  effect  of  angular  spreading,  ¥he  convention  adopted 
for  measuring  ©^  and  ©3  is  that  olookwise  angles  are  positive  and  counter- 
clockwise angles  are  negative. 

The  difference  between  the  percentages  obtained  by  entering  the  upper 
curve  of  Figure  5  with  ©^  and  ©3  is  the  percentage  within  that  range  of 
angles  of  the  spectral  energy  computed  from  the  dispersion  envelope  and 
the  energy  spectrum. 

The  bell-shaped  distribution  of  energy  with  direction.  Figure  5,  was 
inferred  from  limited  evidence  £12"]  and  to  the  knowledge  of  the  authors, 
-  its  exact  form  has  had  limited  observational  verification, 
2G4  Other  Sources  of  Possible  Energy  Loss 

Pierson,  Neumann  and  James  QQ  have  stated  that  in  many  cases, 
dispersion  and  angular  spreading  are  the  only  mechanisms  necessary  to 
account  for  the  energy  decrease  observed  in  wave  spectra  at  locations 
,  remote  from  the  generating  area.  They  held  that  there  may  be  important 


energy  lost  due  to  the  effoot  of  cross-sea  interference,     but  have 
expressed  the  opinion  that  losses  due  to  viscosity  are  generally  negli- 
gible.    The  effects  of  oross   and  opposing  winds  on  wind- gene rated  swell 
arc  not  known,   but  undoubtedly  would  modify  the   spectrum  of  the  swell  to 
an  extont  dependent  on  the  wind  speed  and  the   spectral  composition  of  the 
cross-sea.      It   is  believed  that   at   least  a  qualitative  evaluation  of  these 
effeots  will  bo  possible   as  the  volume  of  test  results  become  greater. 

The  direct  measurement  of  attenuation  presents   serious  difficulties. 
However,   a  careful  analysis  of  the  actual  spectral  energy  at  a  wave 
recorder  as  compared  to  the  forecast   spectral  energy  may  indicate  the    ' 
presence  of  processes  which  cause  energy  losses.     For  instanoe,  theory 
indicates  that  viscosity  effects,    if  they  are   important  at  all,  reduce  the 
energy  primarily  in  the  high  frequency  end  of  the  speotrum  [j3J.      If  high 
frequency  energy  from  a  distant  fetch  is  forecast  to  arrive  at  a  wave 
recorder  at  a  given  time,    and  fails  to  arrive,    it  might  be   inferred  that 
viscosity  effeots  were  responsible  for  the  energy  loss#     This  possibility 
is  discussed  further  in  Chapter  6. 
205     Goals  of  the   Investigation 

The  goal  of  this  paper  is  to  make  available  the  results  of  the 
comparison  of  the  speotral  energy  forecast  by  the  N.  Y.  U.  method  and  the 
spectral  energy  as  obtained  from  a  wave-recording  device. 


1 
The  assumed  effect  of  cross-sea  interference  is  that  the  swell  with 

spectral  components  of  the  same  frequency  as  the  cross-sea  wave  trains 
will  form  unstable   interference  patterns  with  the  latter.     This  will  re- 
sult in  attenuation  of  the   swell  having  speotral  energy  in  those  frequencies. 


It  was  intended  that  these  results  would  permit  conclusions  parti- 
cularly in  the  following  areas : 

1»     The  assumed  form  of  the  energy  spectrum, 

2»     The  assumed  form  of  the  wave-energy     directional  spectrum, 

3*  The  effects  of  viscosity  and  oross-seas  on  the  speotrum  of  the 
swell  waves. 

4#  A  general  verification  of  the  system  whioh  implicitly  includes 
the  items  above. 
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CHAPTER  3 
THE  FORECAST  OF  THE  ENERGY  SPECTRUM 

This  chapter  contains  a  brief  summary  of  the  procedures  followed  by 
the  authors  in  preparing  a  forecast  of  the  energy  spectrum  for  a  given 
observation  time. 

When  a  breeze  comes  up,  the  sea  surfaoe  instantaneously  becomes 
covered  with  tiny  ripples  which  form  more  or  less  regular  arcs  of  long 
radius.  They  increase  rapidly  in  height  as  energy  is  transferred  to  the 
waves  by  the  pushing  and  dragging  forces  of  the  wind.  Soon  a  state  of 
irregular  composite  wave  motion,  called  the  sea,  exists.  The  fully 
developed  state  over  a  long  fetch  with  a  sufficiently  long  duration 
finally  will  be  attained  when  the  wave-generating  traction's  have  increased 
the  total  wave  energy  to  such  a  point  that  dissipation  balances  the  work 
done  by  pushing  and  dragging  forces  on  the  sea  surface.  The  total  amount 
of  wave  energy  accumulated  in  the  composite  wave  motion  is  now  distributed 
over  a  wide  range  of  frequencies.  The  curves  and  assumptions  of  the 
N,  Y,  U,  method  were  used  in  forecasting  the  energy  spectra  in  the  fetch 
area,  H,  0,  603  oontains  several  filters  that  can  be  applied  to  account 
for  the  reduction  in  the  energy  spectrum  as  its  components  move  from  the 
fetch  area  to  an  observation  point.  The  authors,  after  considerable 
investigation  and  study,  oonoluded  that  the  idealized  filters  in  H,  0, 
603  would  not  be  convenient  for  use  in  this  study.  To  achieve  the  de- 
sired results,  a  new  method  of  solving  for  dispersion  and  angular  spread- 
ing was  devised.  The  new  filter  used  when  solving  for  the  effects  of 
dispersion  allowed  for  wind  variations  in  the  fetch  and  resulted  in  a 
continuous  dispersion  envelope  of  periods  present  at  the  forecast  point. 
The  angular  spreading  factor  allowed  for  any  orientation  of  the  fetch, 

11 


3A   Select ion  of  Data 

Appendix  II  describes  the  elaborate  instrumentation  used  by  the 
University  of  California  to  record  continuously  the  water  level  at 
Davenport,  California.  An  examination  of  the  wave  records  available 
showed  that  the  average  wave  heights  were  higher  in  the  winter  at  Daven- 
port due  to  well-developed  storms  in  the  North  Pacific  area.  Therefore, 
the  months  of  January,  February,  and  March  1954  wore  selected  for  closer 
study. 

Weather  maps  at  six-hour  intervals  were  then  studied  for  these  three 
months.  All  possible  wind  circulations  that  resulted  in  fetches  acting 
with  components  toward  Davenport  were  reoorded.  An  envelope  showing  the 
periods  (Ty  and  T^)  and  their  estimated  times  of  arrival  (ETA)  was  com- 
puted for  each  fetch.  These  data  (Ty,  Tl,  ETA)  were  tabled  and  plotted. 
After  detailed  examination,  the  period  between  0630Z,  23  January  1954  and 
0430Z,  1  February  1954  was  seleoted  for  thorough  investigation  and  study. 

This  time  interval  was  the  most  favorable  for  the  following  reasons: 

1*  A  stationary  nearly  uniform  fetch  prevailed.   It  is  discussed  in 
Appendix  III  and  illustrated  in  Figures  7  and  8. 

2.  Representative  wave  reoords  were  available  in  "fast  time"  with  . 
a  time  interval  of  twelve  hours,  as  explained  in  Appendix  II. 

3.  Fetohes  A  and  B  contained  weather  station  ship  "P."  "P"  is 

« 

manned  by  experienced  aerographers.     Therefore,    its  weather  reports  were 
assumed  acourate.     "P"  also  reports   sea-water  temperature    (Ts)  which 

1Figure9/\  is  an  example  of  this  from  0630Z,   23  January  1954  to  0430Z, 
1  February  1954. 

2"Fast  Time"   is  used  to  indicate  the  portions  of  the  record  -where   sea 
level  was  recorded  on  the  scale  of  6  inches  per  1  minute. 

3  See  Appendix  III. 

12 


could  bo  used  in  computing  the  stability  factor,  which  is  used  in 
estimating  the  wind  speed, 

4,  Only  fetches  A,  B,  and  C,  as  defined  in  Appendix  III  and  illus- 
trated in  Figures  7  and  8,  were  believed  to  be  producing  significant 
energy  at  Davenport* 

5,  Looal  winds  in  the  vicinity  of  Davenport  were  blowing  from  the 
land  toward  the  sea  and  were  light. 

6,  All  instruments  were  cleaned  and  calibrated  on  24  January  1954 
and  assumed  to  be  in  good  working  order. 

Considerations  1,  2,   and  4  established  wave  forecast  times  at 
Davenport  which  are  as  follows : . 

tx  043 OZ  28  Jan  1954 

t2  1630Z  28  Jan  1954 

t3  0430Z  29  Jan  1954 

t4  1630Z  29  Jan  1954 

t5  0430Z  30  Jan  1954 

t6  1630Z  30  Jan  1954 

t?  0430Z  31  Jan  1954 

tQ  1630Z  31  Jan  1954 

t9  0430Z  1  Feb  1954 
Hereafter,  the  forecast  time  will  be  referred  to  as  t^, 
3B   Compilation  of  Data 

Table  2  records  the  data  measured  from  the  weather  charts  illustrated 
in  Figure  7# 


•'■See  Table  5. 
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Table  3  records  the  data  oomputed  from  parameters  listed  in  Table  2 
by  the  following  procedures    (symbols  are  all  defined  in  Tables  2  and  3): 

1.  Stability  Factor   (V     N )  depends  on  Ts  minus  T&  and  is  determined 
from  Table  5. 

2.  Isobar  curvature  correction   (l00)  depends  on  isobar  curvature    (I   ) 
and  stability  factor  and  is  determined  from  Table  6« 

3*     Geostrophio  wind  reduced  to  sea  surface   {Yzg) 

Vgs  =   (Vg)    (I00)   (stability  factor)  (4) 

4.     Surface  wind  speed   (V)  obtained  by  oomparing  VQ^,  the  reported 
surfaoe  wind,  with  v*gs* 

5*     Limiting  fetch   (F.  ),  upper  period  CO*   and  E  were  read  directly 
from  Figure  4  using  V« 

6.  Equivalent  duration  (Tg)  was  established  by  comparing  generation 
time  of  last  V  to  present  V,  for  an  equal  E-value« 

7.  Travel  time   (T  t) 


Ro  ^FLG  -  FL 
1.515  T 


T^  r  -^-^ **--  (5) 


Ob    = 


"Where  Ty  leaves  the  wind  area  from  near  the  windward  side  of  the  fetch. 

8.  Lower  period   (T^) 

TL-     ^2 (6) 

L  "       1.515  Tob 

This  Tt  will  have  the   same  ETA  as  Txj  beoause  it  left  from  the  leeward  side 
of  the  fetch. 

9.  Duration  time    (Td) 

Td  ;  TE  -/AT  (7) 

Inhere  /^T  »  6  hours,   or  map  interval. 

See  Figure  7 
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10.     Estimated  time  of  arrival   (ETA) 

ETA  =  map  time  ■/    TOD  (8) 

3C       Development  of  the  Filter 

The  values  of  E*  must  be  reduced  by  filters  as  the  spectral  components 
move  from  the  generating  area  to  the  observation  point*     Chapter  2  explained 
the  theory  of  filters  involved  in  the  N.  Y.  U.  method.     Two  processes  act 
to  limit  the  energy  at  the  observation  point:   dispersion  and  angular  spread- 
ing. 
3D       Dispersion 

Dispersion  results  from  the  presence  of  different  group  velocities 
among  the  spectral  components  in  the  fetch.     The  writers,  to  accurately 
account  for  the  effects  of  dispersion,   plotted  Ttj  and  T*    as  a  function  of 
time. ^     Ty  left  from  the  windward  side  of  the  fetch  with  RQ  />  Fjjj  -  FL 
distance  to  travel*     Tl  moved  out   of  the  leeward  side  with  distance  RQ  to 
travel.     Only  periods  between  Ttj  and  Tt  were  theoretically  present  at 
Davonport  because  of  the  dispersion  influence. 

Fetches  A,   B,   and  C  were  fully-developed  wave  generating  areas ; 
consequently,  upper  period  waves   (Ty)  were  also  leaving  from  the  leeward 
side  with  travel  distance  R0.     Their  ETA  was  computed  by 

T0b  =  "^ 0) 

0D  1.515    Ttj 

and  Equation  (8).     Due  to  a  shorter  travel  distance  this  Ty  arrived  ahead 
of  a  similar  Ty  which  left  from  fetch  windward. 


1See  Table  3 
2See  Table  9A 
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To  include  the  first  arriving  waves  from  the  generating  area,  the 
waves  created  as  the  wind  continued  and  the  waves  leaving  the  fetch  area 
when  the  wind  oeased.  Figure  9A  was  modified  into  Figure  9B,  a  dispersion 
envelope,  A  dispersion  envelope  is  a  plot  of  Ty  and  T^  present  from  the 
same  fetch  as  a  function  of  the  same  ETA.  The  envelope  was  plotted  from 
values  obtained  through  the  use  of  Figure  4  and  Equations  (5),  (6),  (3), 
and  (9).  The  results  are  shown  in  Figure  9B  which  illustrates  the  various 
limiting  envelopes  of  periods  due  to  dispersion  affecting  Davenport  before 
and  during  the  forecast  periods  from  fetohes  A,  B,  and  C* 

The  resulting  reduction  in  E  due  to  dispersion  is  illustrated  in 
Figure  10  for  t^  and  Figure  11  for  t4. 
3E   Angular  Spreading 

Angular  spreading  results  from  the  variability  in  wave  direction  in 
and  after  leaving  the  fetch  areac  The  angular  spreading  factor  is  related 
to  the  orientation  of  the  observation  point  with  respect  to  the  generating 
area  and  is  therefore  dependent  upon  the  parameters  ©3]>  ®4L*  ®3iV*  aia-<^ 
G^jy,  whioh  were  measured  from  the  synoptic  weather  chart.  For  times  near 
the  end  of  a  fetch  life,  important  wave  energy  may  came  from  the  rear  of 
the  fetch  (fotoh  windward).  At  these  times  the  angular  spreading  factor 
decreases  in  value,  consequently,  less  energy  will  be  forecast  to  arrive. 

The  factors  for  the  fetch  leeward  and  fetch  windward  were  obtained 

2 

from  an  angular  spreading  factor  graph.   The  difference  between  these 

two  factors  was  (from  leeward  to  windward)  as  great  as  20,'.  The  authors 
considered  the  possible  error  introduced  by  using  an  angular  spreading 

•'•See  Figure  7 
2See  Figure  5 
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factor  averaged  over  the  entire  fetch  as  unacceptable.  To  accurately 
determine  the  angular  spreading  value  as  a  function  of  distance  from  the 
fetch  front.  Equation  (10)  was  used, 

•  A  D  »  1.515  (Tob)  T  -  R0  (10) 

<&  D  is  the  distance  in  nautical  miles  measured  from  the  leeward  side  of 
the  fetch  into  the  generating  area. 
T0D  ^-s  ^^e  travel  time, 

T  is  any  period  between  Ty  and  Tl  for  t^  as  determined  by  the  dispersion 
envelope, 
R0  is  the  distance  from  observation  point  to  the  center  of  fetch  leeward. 

In  Equation  (10),  T  was  replaced  by  l/frequency,  a  ^D  was  computed 
for  each  0,005  change  in  frequency  for  each  tj_  involving  fetches  A#  B,  and 
C,  The  distance,  ^D  was  thon  calculated  and  plotted,  and  a  corresponding 
G~  and  Q*  measured.  These  angles,  applied  to  the  angular  spreading  graph, 
yielded  the  angular  spreading  factor.  This  procedure  was  usod  when  it  was 
evident  that  the  waves  were  not  leaving  from  the  front  of  the  fetch. 

The  resulting  reduction  in  E  due  to  angular  spreading  is  shown  in 
Figure  10  for  t^  and  Figure  11  for  t^» 

It  is  considered  significant  that  the  use  of  Equation  (10)  allowed 
for  variations  in  angular  spreading  due  to  variations  in  (RQ-/AD). 
Figure  12  is  a  graph  of  angular  spreading  factor  with  respect  to  AD,  for 
fetch  A,  This  factor  did  not  decrease  at  a  constant  rate  in  going  from 
fetch  leeward  to  fetch  windward.  Graphs,  not  included  hero,  yielded 
similar  results  for  fetches  B  and  C,  These  variations  in  rate  of  change 
of  the  angular  spreading  factor  depend  upon  the  orientation  of  the  various 
fetches  with  respect  to  Davenport, 

1-See  Figuro  7 
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3F   Shallow  Water  Effects 

Shallow  water  effects  of  refraction,  shoaling,  percolation,  dif- 
fraction, and  friction  influence  E  as  the  deep  water  wave  energy  approaches 
the  Davenport  area. 

Refraction  was  determined  using  the  method  presented  in  II.  0.  234, 
and  a  wave  refraction  diagram  constructed  using  U.  S.  Coast  and  Geodetic 
Original  Survey  Sheet  No.  5266  of  25  April  1935,  The  chart  shows  the 
bottom  contours  seaward  of  Davenport  to  be  nearly  straight  and  parallel. 
The  range  of  deep  water  wave  directions  was  very  small,  therefore  refraction 
was  treated  as  if  all  deep  water  waves  arrived  from  one  angle.  Shoaling 
effects  were  also  obtained  from  H.  0.  234,  Plate  I. 

Percolation  and  friction  were  assumed  negligible  due  to  the  narrow 
width  of  the  continental  shelf.  Diffraction  was  assumed  to  be  negligible. 
All  instruments  were  located  in  open  water  well  away  from  rock  formations 
or  man-made  structures. 

The  shallow-water  effeots  of  refraction  and  shoaling  wore  combined 
into  one  correction  factor.  Their  influence  on  E  is  shown  in  Figure  10 
for  t^  and  Figure  11  for  t4»  Figures  10  and  11  show  the  resultant  ocean 
wavo  spectral  energy  forecast  to  arrive  at  the  recording  instruments  for 
t^  and  t^..  The  area  under  the  curve  can  be  integrated  by  a  planimeter 
to  evaluate  the  total  energy  present  for  any  forecast  time.  The  results 
aro  shown  in  ohapter  5,  and  are  plotted  with  respeot  to  time  in 
Figure  14. 
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CHAPTER  4 
SPECTRUM  ANALYSIS 

The  objective  of  this  chapter  is  to  discuss  the   general  procedure 
and  theory  employed  in  obtaining  the  energy  spectrum  from  the  record  of  a 
wavo-reoording  device*  ; 

In  this   study,  wave  reoords  were  used  from  both  a  surface  wave  record- 
ing device  and  a  pressure  recording  device.       The  distinction  between  the 
two  is  that  the  former  reoords   surface  waves  and  the  latter  records  wave 
pressures  at  the  depth  in  which  it  is   installed* 
4A       Theory 

Many  of  the  features  of  surface,    and  below  surface  waves  have  been 
explained  by  Pier  son  and  Marks   OQ  and  Neumann   [loj  based  on  the  assumption 
that  a  wave  record  as  a  function  of  time  is  the  rosult  of  a  Gaussian 
process.     The  theory  was  developed  by  Tukey  and  Hamming  [_9]  among  others. 
In  the  publication  cited,  the  latter  explain  a  method  which  may  be  used 
for  estimating  the  power  spectrum  of  a  steady-state  sea  surface  using  the 
Fourier  transform  of  the  autocorrelation  function.     A  text  which  treats 
correlation  functions  and  power  spectra  in  detail  is   Information  Theorynz] 
by  Goldman. 

Tirame  and  Stinson   [l4J  present  curves  of  the  pressure  response 
factor  as  a  function  of  wave  period,  T,   and  depth,  H.     The  pressure 
response  .factor   is  a  measure  of  the  attenuation  of  wave  energy  with  depth. 
4A1     numerical  Estimation  of  Non-normalized  Autocorrelation  Function      1 

and  the  Pressure  Powor  Spectrum 

The  results  of  Tukey  and  Hamming  £sQ  are  used  to  derive  the  power 

-"■See  Appendix  III 
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spectrum  from  a  pressure  record  placed  at  the  sea  floor  in  a  depth  of 
water,  II.  The  possible  errors  which  enter  into  the  numerical  analysis 
have  been  considered*  The  corrections  which  must  be  applied  to  the  raw- 
results  and  the  degree  of  reliability  have  been  determined. 

The  formulas  for  numerical  estimation  of  the  power  spectra  as 
presented  by  Tukey  and  Hamming Ql  and  as  used  by  Pierson  and  Harks  {8}    are 
given  by: 

*Z   "  ^l*  =  t3  "  "^2  = %  "  %-l  a  *£&  (11) 

N-P 

%  '   iTTp  />        P(tn^  P^n-P)   0*0.1.2,  —  m)     (12) 
n=l 

m-1 

L    :-    (Qo/2        >  Qp  Gos7^^-     -/ Qfa.  CosTTh)      (13) 

m  '      »  m 

P=l 

(h»0.1, m) 

Uh  -  .23Ih-1/.54  lfc/,23  hki^1   (h*0,lf m)  (14) 

[A^)]2.'-^  (16) 

Where/<  -  2  TT  f 

Equation  (11 )  states  that  N-l  equally  spaced  time   intervals  are 
marked  off  on  the  record  by  points  spaced  £*t  units  apart*     The  values 
of  pressure  at  these  points  in  terms  of  departures  from  the  mean  are  then 
tabulated  as  P(ti),   P(t2) P("tn), p0*n).  h  numbers  resulting. 

The  non-normalized  lag  coefficients  are  then  found  by  the  use  of 
Equation   (12),      Qp  ranges  from  Qq  to  Qg^   or  m  •/  1  total  numbors.     Qq  is 
the  Ep^  max,   or  the  total  spectral  energy  uncorrected  for  the  scale  faotor 
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of  the  record,   at  the  depth  II. 

The   summation  of  Equation  (12)  was  accomplished  by  the  use  of  a 

CRC-102-A  computer.1 

The  "raw"  estimates  of  power  in  the/A  band  from/d  =        \"~?J  to 

Z>  tm 

^  miiY^rS)     are  found  from  Equation  (13).  A  plot  of  the  "raw"  estimates 
with  values  of  h  result  in  a  quite  jagged  graph.  The  jaggedness  is  intro- 
duced by  the  measurement  of  P  at  disorete  points  and  by  the  finite  length 
of  the  record.  A  smoothing  device  is  given  by  Equation  (14). 

The  final  numbor,  Uh#  is  an  estimate  of  the  mean  value  of  the  power 

system  from  "ft1  "  ?•<     to  *^-%rt~  •  A"b  "bhe  point  given  by  Equation  (15) 
^\tm  *•*  tan 

the  best  estimate  of  the  value  of  IK^  (^fo)J  *-s  Sivon  by  Equation  (16). 
The  result  is  a  smoothing  of  the  "jagged"  curve  of  Equation  (13),  and  this 
continuous  curve  is  the  best  estimate  of  the  pressure  power  spectrum. 

For  the  purposes  of  this  investigation,  the  computer  program  for  the 
Fourior  analysis  was  limited  to  Equation  (13).  The  refinements.  Equations 
(14),  (15),  (16)  were  not  included  in  the  solution,  therefore,  the  resulting 
spectra  are  the  "raw"  or  unsmoothed  curves.  This  limitation  was  dictated 
by  time  considerations  as  the  preparation  of  such  a  program  is  extremely 
lengthy.  These  refinements  will  be  included  in  the  computer  program  as 
time  permits  and  will  be  available  for  similar  investigations  in  the 
future. 
4A2  Reliability  of  Results 

The  final  results  of  paragraph  4A1  are  estimates*  Estimates  can  be 
in  error,  and  the  important  part  of  this  analytical  method  is  that  it 
also  presents  an  estimate  of  the  error.  Tukey  and  Hamming  [_9j  have  shown 

jDhe  computer  programming  of  the  equations  used  in  this  investigation 
was  accomplished  by  Professors  J.  B.  Wickham  and  H.  M.  Martinez  of  the 
U.  S.  IJaval  Postgraduate  School. 
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'  2 

that  the  errors  follow  a^C  distribution  with  f  degrees  of  freedom  given 


f  s  L:l*  (17) 


by 


Givon  f,  thon  the  use  of  Table  (7)  yields  information  on  the  reliability 
of  the  estimate  for  each  frequency  band  obtained  by  the  analysis. 

For  20  degrees  of  freedom,  theUv  obtained  Cram  Equation  (16)  may  be 
multiplied  by  1*8  and  by,63.  The  true  value  of  U  in  the  band  under  analysis 
will  lie  between  l,8Uh  and  ,63U  ^  90^  of  the  time, 
4B   Procedure 

For  this  investigation,  the  verification  periods  were  of  the  order  of 
twenty  minutes  on  the  wave  record.  The  &t   was  chosen  as  two  seconds  as 
determined  by  the  criterion  of  Pierson  and  Harks  jjQ  •  The  total  number 
of  measurements  per  wave  record,  N,  was  of  the  order  of  GOO*  For  pro- 
gramming convenience,  the  total  number  of  lags,  m,  was  taken  as  43, 

The  total  number  of  consecutive  verifications,  12  hours  apart,  was 
nine.  It  was  originally  intended  that  an  autocorrelation  program  and  a 
Fourier  analysis  would  be  carried  out  for  each  of  these  times.  However, 
some  initial  difficulties  in  the  programming  routine  restricted  the  time 
available  for  the  investigation,  and  time  did  not  permit  the  fulfillment 
of  this  objective. 


In  future  work,  the  wave  record  should  be  analyzed  over  a  longer  time 
interval  resulting  in  greater  value  of  IT.  This  will  increase  the  degrees 
of  freedom  and  permit  better  resolution  of  the  spectral  curve.  The 
number  of  lags  taken  should  be  increased, 

2 
See  pp  13,  Chapter  3. 
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CHAPTER  5 
STATEMENT  OF  RESULTS 

The  purpose  of  this  chapter  is  to  present  the  results  of  this 
investigation, 
5A   The  "Filter" 

It  is  the  opinion  of  the  authors  that  the  "filter"  devised  for  the 
energy  spectra  forecast  in  this  study  represents  an  accurate,  practical 
method  of  solution.  This  is  particularly  true  for  the  case  of  marked 
variation  of  wind  speed  with  time  in  a  generating  area.  For  example,  in 
the  case  of  increasing  wind  speed  with  time,  energy  fronts  generated  by 
the  earlier,  lower  speed  winds  will  frequently  be  overtaken  by  energy 
fronts  generated  by  the  later,  higher  speed  winds.  This  additional 
energy,  from  the  higher  speed  wind  might  be  overlooked  if  a  plot  of  the 
"dispersion  envelope"  with  time  were  not  maintained. 

The  filters  and  examples  described  in  the  N.  Y.  U.  method  6 
primarily  treat  the  simplified  cases  of  constant  or  slowly  changing  wind 
speed.  This  permits  the  energy  to  be  obtained  directly  from  the  G.G.S. 
curves  for  one  wind  speed* 

It  must  be  understood  that  the  "period  envelope"  must  be  identified 
precisely  with  the  wind  speed  -which  generated  it.  It  appears  that  the 
filter  herein  described2  adequately  handles  the  more  complicated 
situations. 

It  appears  that  in  forecasting  by  this  method,  using  the  filter 
presented  herein,  the  forecaster  will  frequently  need  to  work  directly 
with  the  spectral  curves  in  addition  to  the  C.C.S.  curves.  For  this 


See  Figure  9B 


2See  Chapter  3 
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reason,  the  tabled  solution  to  Equation  (1)  is  included  in  Appendix  I. 
It  is  recommended  that  personnel  using  this  method  plot  the  curves  of 
the  spectra  for  various  wind  speeds  on  an  enlarged  scale  for  the  purpose 
of  representing  the  spectrum  present  at  the  locality  and  time  of  interest, 
5B  Verification  of  the  Forecast 

Spectrum  forecasts  were  prepared  for  Davenport,  California,  for  the 
times  ti  through  tq  as  described  in  Ckapter  3.  These  forecasts  were 
prepared  for  two  geographical  sites  having  water  depths  of  50  and  81.5 
feet  respectively.  A  description  of  these  sites  and  the  instrumentation 
of  each  is  presented  in  Appendix  II. 

Data  available  for  verification  of  these  forecasts  were  the 
following: 

1.  Significant  wave  height  measurements  for  the  sea  surface-water 

depth  50  feet,  hereafter  called  site  I.  These  data  were  computed  from 

sea  surfaoe  wave  records  by  personnel  of  the  University  of  California, 

Berkeley,  and  made  available  to  the  authors  for  this  study.  Significant 

wave  height,  H  i  ,  is  defined  as  the  average  of  the  heights  of  the  one- 
s' 
third  highest  waves.  E  values  may  be  reoovered  by  the  statistical 

relationship, 

-  =  2.83^  (18) 

o 

A  plot  of  the  E  values  obtained  in  this  manner  for  the  times  t^  through 
to,  is  illustrated  in  Figure  14. 

2.  A  sea-level  wave  record  for  site   I  obtained  from  a  Beach  Erosion 
Board,   Step-resistor  gauge  as  described  in  Appendix  II.     These  records 
were  analyzed  by  the  procedures  described  in  Chapter  4,   and  an  estimate 

1See  Figure  13. 
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of  the  E  value  for  times  ti  and  t2  were  obtained. 

3.     A  "pressure"  reoord  for  the  site  in  the  water  depth  of  81.5  feet. 
This  site  will  be  referred  to  as   site   II.     This  pressure  reoord  was  ob- 
tained from  a  Mark  IX  pressure-head  instrument  described  in  Appendix  II. 
The  instrument  was   installed  1%  feet  from  the  bottom  at   site   II  or  in  a 
water  depth  of  approximately  80  feet.     These  records  were  analyzed  by  the 
procedures  described  in  chapter  4  and  an  estimate  of  the  E-value  for 
times  t3  through  tg  were  obtained. 

As  described  above,  the  "records"  available  were  obtained  from  two 
different   instruments.     For  the  times  t^  and  tg,  the   surface   step- 
resistance  devioe  was  available  at   site   I.     For  the  times  t3  through  tg 
the  pressure-head  device,   installed  at  a  depth  of  80  feet,  was  available 
at   site   II. 

Figure  14  is  a  plot  of  E- values  for  the  times  t^  through  tg  for  the 
sea  surface  at   site   I,   as  estimated  by; 

1.  The  forecast. 

2.  The  H  i  computations  provided  by  the  University  of  California, 

3" 
as  described  above.  The  90^  confidence  limits  IjQ  are  plotted  at  each 


point. 

3.     The  writers'   analysis  of  the  water  level  record  of  the   step- 
resistance  gauge  as  described  above.     The  90^0  confidence  limits  for 
these  points  are  also  included. 

Figure  15  is  a  plot  of  E-value s  for  the  times  tg  through  tg  for 
site   II  for  a  depth  of  80  feet,   as  estimated  by: 

1*     The  forecast. 

^•See  pp  22,   Chapter  4. 
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2.  The  -writers'  analysis  of  the  "pressure"  record  of  the  instrument 
as  described  above.  The  90^  confidence  limits  for  these  points  are  also 
included,-*- 

The  prooedure  for  obtaining  the  foreoast  spectra  and  E-values  for 
verification  times  t3  through  tg  at  site  II  was  the  following, 
A  surfaoe  spectrum  forecast  for  a  given  verification  was  prepared  at  site 
II,  This  spectrum  was  reduced  to  a  depth  of  80  feet  using  the  pressure 
response  attenuation  factors  of  Timme  and  Stinson  Q.<T]  as  described  in 
Ghapter  4,  An  example  of  the  results  of  such  a  prooedure  is  illustrated 
in  Figures  16  and  17  for  verification  times  tg  and  ty  respectively. 

As  stated  earlier,  the  original  intention  was  to  obtain  the  computed 
spectrum  for  each  verification  time.  However,  because  of  some  initial 
difficulties  onoountered  in  the  computer  programs  discussed  in  Ghapter  4, 
the  time  allotted  for  the  use  of  the  computer  did  not  permit  the  fulfill- 
ment of  this  objective, 

A  disoussion  of  the  figures  described  above,  and  the  conclusions  of 
the  authors  will  be  presented  in  Ghapter  6, 

*See  pp  22,  Ghapter  4. 
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CHAPTER  6 
CONCLUSIONS 

The  purpose  of  this  chapter  is  to  oompare  the  observations  with  the 
forecasts  and  to  suggest  possible  explanations  for  any  differences  noted. 
6A   General  Comments  on  the  Results 

An  examination  of  Figures  14  and  15,  ■which  compare  forecast  and  ob- 
served energy,  and  Figure  22,  which  compares  the  forecast  and  observed 
significant  height,  indioates  certain  interesting  features. 

1*   It  is  noted  that  the  general  agreement  is  quite  good.  The  fore- 
cast energy  peaks  oocurred  at  the  time  predicted  and  the  energy  was  of  the 
correct  order  of  magnitude* 

2,  The  forecast  energy  was  in  general  equal  to  or  greater  than  the 
observed, 

3*  At  the  point  of  greatest  difference  the  forecast  significant 
height  was  about  one  and  one-half  feet  greater  than  the  observed.  This 
would  indicate  that  for  these  conditions,  the  forecast  would  be  of 
operational  value, 

4#  It  is  noted  that  the  points  of  closest  agreement  occur  where  the 

observed  energy  is  at  or  near  a  peak, 

6B   Possible  Explanations  for  the  Differences  between  the  Observed  and 
Forecast  Energy 

For  the  purpose  of  suggesting  possible  explanations  of  differences 
between  the  observed  and  forecast  energy,  the  following  items  are  con- 
sidered to  be  possible  sources  of  error: 

1,  The  theory.  The  foreoast  energy  depends  upon 

a,  the  energy  spectrum  and, 

b,  the  directional  spectrum. 
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If  the  assumed  form  of  either  of  those  curves  is  incorrect,  the  energy 
forecast  would  probably  be  in  error. 

2.  The  application.   In  applying  the  method,  it  is  unavoidable  that 
the  selection  of  data  is  largely  subjective.   It  is  possible  therefore  that 
sizable  errors  could  be  introduced.  Errors  in  the  following  mi[;ht  be 
particularly  critical: 

a*  geostrophic  wind  analysis 

b#  surface  wind  calculation 

c.  measurement  of  fetch  parameters. 

3.  Energy  transformation  processes.  Processos  or  conditions  which 
are  believed  to  affect  ocean  swell  wave  energy  and  are  not  accounted  for  in 
the  forecast  are  the  following: 

a.  oross-sea  interference 

b.  cross  or  opposing  wind 

c.  viscosity  effects. 

4.  Refraction  and  shoaling  effects. 
6C   The  Energy  Forecast 

Figure  14  oompares  the  forecast  value  of  E  with  the  observed  at  site  I 
at  the  sea  surfaoe  for  the  verification  times  t^  through  tg.  It  is  noted 
that  the  foreoast  E  value  lies  above  the  90^o  confidence  limits  of  the 
observed  E  value  at  the  times  t,,  t,,  tg,  and  tg.  The  inference  might  be 
drawn  that  at  these  times  energy  loss  processes,  not  considered  in  the  fore- 
oast, are  of  importance*  It  is  believed  significant  that  a  well  defined 


Additional  computations  of  these  data  indicate  that  the  correct  value  of 
E  is  approximately  identical  with  the  foreoast  value. 
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aroa  of  oross-soa  was  present  at   such  location  and  times  as  to  have  inter- 
fered with  the   "swell"  arriving  at  the  times  tg,  t^,   and  tg  from  fetohos  A 
and  B   (Appendix   III).     Gross-seas  are  believed  to  oauso  the   summation  of  the 
spectral  components  to  become  non-linear  giving  rise  to  the  conversion  of 
wave  energy  to  turbulent  onorgy  by  formation  of  white-caps.     Figure  18 
illustrates  the  looation  and  orientation  of  the  cross-sea  rolativo  to  the 
generating  fetohes  and  the  foreoast   sites  at  Davenport.      Figure  19  displays 
the   spectrum  forecast  for  the  cross-soa  generating  area  at   six  hour  inter- 
vals  over  the  period  when  it  is  most  clearly  defined.     For  reasons  previously 
stated,  oomputed  spectra  for  the  verification  times  under  discussion,   are 
not  available.      It   is  believed  that   such  spectra  would  have  aided  greatly  in 
the  determination  of  the  portion  of  the  spectrum  affected  by  the  cross-sea 
interference.     However,    it  was  noted  that  if  the  foreoast   swell  energy  were 
tobally  erased  in  the  frequency  band  of  the  cross-sea3,  energy  reductions 
of  the  order  of  those  observed  in  Figure  14  would  have   occurred. 

Figure  15  oomparos  the  forecast  values  of  E  with  the  observed  values 
at  site   II,   at  a  depth  of  80  feet,   for  the  times  t3  through  tg. 

By  comparing  the  energy  differences  which  ocour  at  the  surface 
(Figure  14)  with  the  energy  differences  whioh  ooour  at  depth   (Figure  15) 
it  appears   it  might  be  permissible  to  draw  inferences  about  the  frequency 
ranges  in  whioh  these  differences  occur.     These   inferences  are  based  on 
the  fact  that  the  middle  and  lower  frequency  energy  would  be  less 
attenuated  with  depth  by  the  hydro-dynamio   filter  than  would  the  high 
frequency  energy.     Accordingly,   if  the  relative  error  of  the  energy  differ- 
ences at  the  surface  was  of  the  same  order  as  the  relative  error  at  depth, 
it  might  be  inferred  that  the  energy  differences  wore  primarily  in  the 
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middle  and  lew  frequencies.  By  similar  reasoning,  if  the  relative  error 
at  depth  was  significantly  less  than  the  relative  error  at  the  surfaoe, 
it  might  be  inferred  that  the  energy  differences  were  primarily  in  tho  high 
froquenoy  end  of  the  spectrum.   If  the  relative  error  at  depth  is  greater 
than  the  relative  error  at  the  sea  surface  the  possible  causes  are  less 
obvious.  A  possible  explanation  might  be  that  too  much  middle  or  low 
frequency  energy  and  too  little  high  frequency  energy  is  forecast,  the 
errors  cancelling  and  giving  close  agreement  with  the  observed  energy  at 
the  surface.  Then,  as  tho  forecast  spectrum  is  reduced  by  the  hydro- 
dynamic  filtor,  at  depth  tho  surplus  middle  and  low  frequency  energy  which 
was  forecast  would  be  relatively  less  attenuated  with  depth  than  the  high 
frequency  energy.  Then  at  depth,  the  forecast  energy  by  virtue  of  the 
surplus  middle  frequency  energy  would  show  a  greater  rolative  error  than 
at  the  surfaoe. 

It  should  be  noted  that  it  is  quite  possible  that  minor  wind  systems 
suoh  as  sea-breeze  circulations,  may  have  developed  between  the  six  hour 
map  intervals  of  this  study  which  could  have  developed  energy  sources  not 
taken  into  account.  Purthor,  beoause  site  II  is  some  distance  seaward  of 
site  I,  tho  total  energy  at  sito  II  is  greater  as  tho  shoaling  effect  is 
less  pronounced.  Since,  in  the  above  discussion,  the  energy  differences 
are  compared  between  site  I  at  the  sea  surface  and  site  II  at  depth,  the 
relative  error  as  observed  from  Figure  15  will  be  somewhat  greater  than 
would  be  the  rolative  error  at  depth  at  site  I. 

Examining  the  relative  errors  of  energy  differences  at  the  sea- 
surface  and  at  depth  for  the  times  tj  through  tn,  and  applying  the 

Refraction  differences  betweon  site  I  and  site  II  wore  assumed 
negligible. 
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roasoning  suggested  above,  the  following  comments  are  offered: 

1*  At  the  times  t%   and  t5  the  relative  erros  at;  the  sea-surf aoo  and 
at  depth  are  about  the  same.  At  the  time  t^  the  relative  error  is  much  less 
at  depth  than  at  the  surface,  Following  the  roasoning  suggested  above  it 
might  be  inferred  that  at  times  tg  and  tg  important  energy  differences 
occurred  in  the  middle  frequencies  and  that  at  time  t^  the  energy  difference 
was  primarily  in  the  high  frequency  range.  As  noted  previously,  it  is 
believed  that  oross-sea  interference  acted  on  the  "swell"  arriving  at  these 
times*  These  comparisons  give  rise  to  the  possible  conclusion,  that  the 
foreoast  spectrum  ovorprediots  in  the  low  and  middle  frequencies  and  under- 
predicts  in  the  high  frequencies," 

2.  At  the  times  tg  and  ty,  the  relative  error  of  the  energy  differ- 
ences was  significantly  greater  at  depth  than  at  the  surface.  As  suggested 
above,  a  possible  explanation  for  this  may  be  an  overforecast  of  enorgy  in 
the  middle  and  low  frequencies  and  an  underforeoast  of  energy  in  the  high 
frequency  band.   It  may  be  noted  that  these  times  are  at  or  near  an  energy 
peak.   It  might  be  suggested  than  that  for  higher  wind  speeds  the  assumed 
enorgy  spectrum  overprodiots  in  the  middle  frequencies  and  underpredicts 

in  the  higher  frequencies.  At  least  suoh  an  assumption  is  consistent  with 
these  data. 

3.  At  the  times  to  and  to,  the  relative  error  at  the  surface  is  signi- 
ficantly larger  algebraically  than  at  depth.  Following  the  line  of  reason- 
ing used  for  time  t^,  it  might  be  inferred  that  important  energy  differences 
occurred  in  the  high  frequoncy  band.  The  oause  of  this  assumed  high 
frequency  energy  reduction  is  not  readily  apparent.   It  is  noted,  however, 
that  at  both  time  t^.  and  tg  the  observed  energy  is  at  a  minimum.  A 
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possible  explanation  might  be  that   at   lower  wind  speeds  the   assumed  enorgy 
spocbrum  may  overprediot   in  the  high  frequency  end  of  the  spectrum. 
6D       Summary 

In  summary,    it   is  the  opinion  of  the  authors  that  the  forecast  method 
for  swell  is  operationally  useful  for  fotohes  of  the  type  active   in  this 
study.     Lased  on  this   investigation,   it  appears  that  the  ourves  and  pro- 
cedures of  the  method  fjQ  ,  when  properly  applied,  will  result  in  a  wave 
forecast  of  useful  accuracy.     Figure  22,   a  plot  of  the  forecast   significant 
wave  height  compared  to  the  observed  illustrates  this  conclusion*     As  noted 
in  paragraph  6A  above,  the  results  are  consistent  with  the  assumption  that 
energy  reduction  processes  not  considered  in  the  forecast  are  at  times 
important.      In  particular,    it   is  believed  that  cross-sea  interference  was 
responsible  for  significant  energy  reductions   at  certain  verification  tunes 
in  this   investigation.      It   is  believed  that  a  qualitative  correction  would 
be  possible  by  reducing  tho  forecast   spectrum  by  a  fraction  of  the  cross-sea 
speotrum  developed. 

The  data  are  also  consistent  with  tho  premise  that,  for  increasing 
wind  speeds,  the  forecast  energy  speotrum  is  too  high  in  the  low  and  middle 
frequencies  and  possibly  too  low  in  the  high  frequency  band#   and  that  for 
lower  wind  speeds,   an  opposite  trend  appears  possibleo 
6E       Future  Research 

When  ce**tain  modifications   in  the  computer  programs  discussed  in 
Chapter  4  are  completod,   and  accurate  observational  spectra  are  available, 
more  dofinito  conclusions  will  be  possible.      It   is  believed  that  cross-sea 
effects  may  beoome   quite  clear  from  a  comparison  of  the  forecast  and  the 
observed  spectra.      It   is  the  opinion  of  the  authors  that  extensions  of  the 
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methods  usod  in  this  investigation  can  provide  the  quality  and  quantity 
of  data  nooessary  for  the  evaluation  of  forecast  methods  and  the  effects 
of  attenuation. 
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APPENDIX  I 
The  purpose  of  this  appendix  is  to  tabulate  values  of  fk   (f  )y-   in 
Equation  (l).  Chapter  1,  as  a  function  of  wind  speed  and  frequency. 


.2 


£(f^2 -TT^Ts  e  "fw  (1) 

(2  7T)&f°  4TT  f  v 

where  C  -  4«8  x  10^  cm^  sec~5 

The  values  of  the   /A(f)l      for  the   solution  of  Equation  (la)  are  given  in 
Table    (8). 
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APPENDIX  II 

DESCRIPTION  OF  LOCALE  AND  FACILITIES  AT 
THE  FORECAST   SITE,    DAVENPORT,    CALIFORNIA 

The  University  of  California,  with  the  permission  of  Santa  Cruz 
Portland  Cement  Company,   planned,   established,   and  maintained  an  ocean 
wave  recording  program  at   Davenport,   California,     As  shown  in  Figure  15, 
the  pier  extended  twenty-two  hundred  feet   seaward  of  the  beach  and  was 
exposed  to  the  large  waves  from  the  North  Pacific  generating  areas. 
Elaborate  instrumentation  and  recording  devices  wero  used.     Continuous 
ocean  wave  records  from  8  Nov  1952  to  1  Dec   1954  were  obtained  from  the 
Beach  Erosion  Board   (BEB),   step-resistor  gages    (parallel  type)  £^15}  and 
Mark  IX  [1Q  pressure  heads.     The  locations  of  these  two  instruments  are 
shown  in  Figure  13.     The  BEB  was  looated  at  the   seaward  end  of  the  pier. 
Depth  immediately  below  the  BEB  was  50  feet.     The  Mark  LX  was  located  in 
8li  feet  of  water   (below  MLLW). 

Fluctuations   of  the  sea  surface  were  transferred  by  electrical  cir- 
ouitry  to  a  recording  center  located  on  shore.     Here  the  fluctuations  were 
permanently  recorded  on  Esterlino-Angus  Continuous  Record  Charts.       Two 
time   scales  were  used: 

6"  per  1  minute,   called  "fast  time" 
6"  per  1  hour,   called  "slow  time" 
Fast  time  was  recorded  every  twelve  hours  for  twenty  minutes  beginning  at 
0430Z  and  1630Z  of  eaoh  day.     tj_   (forecast  time)v*as  determined  by  the  fast 
time  records. 

A  continual  maintenance   schedule  was  kept,  which  involved  mainly 
cleaning  and  re-calibration. 


1See  Figure  21 
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APPENDIX  III 

DESCRIPTION  0?  THE  SYNOPTIC  SITUATION 
DURING  AND  PRECEDING  THE  FORECAST  TEE 

As  explained  in  Chapter  3,  tho  forecast  time  ber;an  at  0430Z,  28  Jan 
1954  and  ended  at  0430Z,  1  Feb  1954. 

The  synoptio  situation  affecting  the  forecast  locality  during  and 
preceding  the  forecast  time  remained  basically  stationary.  A  low  pressure 
center  (average  1000  MB)  existed  over  the  Gulf  of  Alaska  or  British 
Columbia  from  1230Z,  22  Jan  1954  to  1230Z,  28  Jan  1954.  A  strong  high 
pressure  center  (1040  MB)  maintained  itself  over  Alaska  and  the  Bering  Sea 
with  a  ridge  extending  southeastward  down  to  35  N  latitude.  Figure  7 
illustrates  the  general  synoptic  situation  that  existed  before  and  remained 
nearly  constant  during  the  forecast  period. 

The  large  cyclonic  circulation  in  the  Gulf  of  Alaska  created  a  fetch 
of  length  (Flq)  1320  nautical  miles.  Fetch  width  (F^)  was  520  nautical 
miles  v^ith  the  distance  from  the  center  of  tho  leeward  side  to  Davenport 
(R  )  being  420  nautical  miles.  This  fetch  hereaftor  will  be  referred  to 
as  Fetch  "A."  Fetch  A  existed  fully  developed  from  0030Z,  23  Jan  1954  to 
0630Z,  26  Jan  1954. 

At  approximately  0630Z,  26  Jan  1954,  due  to  a  slight  reorientation  of 
the  cyclonic  circulation,  Fj^q.  became  800  nautical  miles  with  a  correspond- 
ing inorease  in  R0  to  960  nautical  miles.  This  change  resulted  in  a 
redesignation  to  Fetch  "B."  All  other  geographic  variables  remained 
constant.  Fetch  B  existed  from  0630Z,  26  Jan  1954  to  1230Z,  28  Jan  1954. 

Data  for  fetches  A  and  B  were  recorded  in  Table  2,  Chapter  3,  and 
the  fetches  are  illustrated  in  Figure  7. 

Wind  is  one  of  the  most  important  variables  in  forecasting  an  energy 
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spectrum.  The  general  wind  area  remained  fairly  constant  bub  wind  varied 
before  and  during  the  forecast  period  in  fetches  A  and  B*  The  variation 
of  the  average  wind  speed  with  time  in  fetches  A  and  B  is  illustrated  in 
Figure  17* 

Figure  8  illustrates  the  general  synoptic  situation  that  existed  from 
0030Z,  27  Jan  1954  to  1830Z,  27  Jan  1954.  A  fetch,  Fetoh  C,  was  fitted 
to  this  circulation.  Data  resulting  from  this  fetoh  are  recorded  in 
Table  2,  Chapter  3, 

The  wind  field  remained  constant  at  22  knots  during  the  existence  of 
Fetoh  C. 

Fetch  A  existed  for  84  hours,  B  for  48  hours,  and  C  for  18  hours.  The 
subscript  notation  (A.,B^,Cj )  was'used  to  distinguish  fetch  sequence,  with 
a  map  interval  of  6  hours.  The  fetch  life  started  when  i  =  1,  and  ended 
when  i  »  14,  i  »  8,  and  i  »  4  respectively. 
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